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(57) Abstract: Techniques for processing a data transmission at the transmitter and receiver. In an aspect, a time-domain implemen- 
tation is provided which uses frequency-domain singular value decomposition and "water-pouring" results to derive time-domain 
pulse-shaping and beam-steering solutions at the transmitter and receiver. The singular value decomposition is performed at the 
transmitter to determine eigen- modes (i.e., spatial subchannels) of the MIMO channel and to derive a first set of steering vectors 
used to "precondition" modulation symbols. The singular value decomposition is also performed at the receiver to derive a second set 
of steering vectors used to precondition the received signals such that orthogonal symbol streams are recovered at the receiver, which 
can simplify the receiver processing. Water- pouring analysis is used to more optimally allocate the total available transmit power to 
the eigen-modes, which then determines the data rate and the coding and modulation scheme to be used for each eigen-mode. 
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TIME-DOMAIN TRANSMIT AND RECEIVE 
PROCESSING WITH CHANNEL EIGEN-MODE 
DECOMPOSITION FOR MIMO SYSTEMS 

BACKGROUND 

Field 

[1001] The present invention relates generally to data communication, and more 
specifically to techniques for time-domain transmit and receive processing with channel 
eigen-mode decomposition for multiple-input multiple-output (MIMO) communication 
systems. 

Background 

[1002] In a wireless communication system, an RF modulated signal from a 
transmitter may reach a receiver via a number of propagation paths. The characteristics 
of the propagation paths typically vary over time due to a number of factors such as 
fading and multipath. To provide diversity against deleterious path effects and improve 
performance, multiple transmit and receive antennas may be used. If the propagation 
paths between the transmit and receive antennas are linearly independent (i.e., a 
transmission on one path is not formed as a linear combination of the transmissions on 
other paths), which is generally true to at least an extent, then the likelihood of correctly 
receiving a data transmission increases as the number of antennas increases. Generally, 
diversity increases and performance improves as the number of transmit and receive 
antennas increases. 

[1003] A multiple-input multiple-output (MIMO) communication system employs 
multiple (Nt) transmit antennas and multiple (N R ) receive antennas for data 
transmission. A MIMO channel formed by the Nr transmit and Nr receive antennas 
may be decomposed into Nc independent channels, with Nc ^ min {Nt, Nr) . Each of 
the Nc independent channels is also referred to as a spatial subchannel of the MIMO 
channel and corresponds to a dimension. The MIMO system can provide improved 
•performance (e.g., increased transmission capacity) if the additional dimensionalities 
created by the multiple transmit and receive antennas are utilized. 
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[1004] The spatial subchannels of a wideband MIMO system may experience 
different channel conditions (e.g., different fading and multipath effects) across its 
bandwidth and may achieve different signal-to-noise-and-interference ratios (SNRs) at 
different frequencies (i.e., different frequency bins or subbands) of the overall system 
bandwidth. Consequently, the number of information bits per modulation symbol (i.e., 
the data rate) that may be transmitted at different frequency bins of each spatial 
subchannel for a particular level of performance may be different from bin to bin. 
Moreover, the channel conditions typically vary with time. As a result, the supported 
data rates for the bins of the spatial subchannels also vary with time. 
[1005] To combat the frequency selective nature of the wideband channel (i.e., 
different channel gains for different bins), orthogonal frequency division multiplexing 
(OFDM) may be used to effectively partition the system bandwidth into a number of 
(N F ) subbands (which may be referred to as frequency bins or subchannels). In OFDM, 
each frequency subchannel is associated with a respective subcarrier upon which data 
may be modulated, and thus may also be viewed as an independent transmission 
channel. 

[1006] A key challenge in a coded communication system is the selection of the 
appropriate data rates and coding and modulation schemes to be used for a data 
transmission based on channel conditions. The goal of this selection process is to 
maximize throughput while meeting quality objectives, which may be quantified by a 
particular frame error rate (FER), certain latency criteria, and so on. 
[1007] One straightforward technique for selecting data rates and coding and 
modulation schemes is to "bit load" each frequency bin of each spatial subchannel 
according to its transmission capability, which may be quantified by the bin's short- 
term average SNR. However, this technique has several major drawbacks. First, 
coding and modulating individually for each bin of each spatial subchannel can 
significantly increase the complexity of the processing at both the transmitter and 
receiver. Second, coding individually for each bin may greatly increase coding and 
decoding delay. And third, a high feedback rate may be needed to send channel state 
information (CSI) indicative of the channel conditions (e.g., the gain, phase and SNR) 
of each bin. 
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[1008] There is therefore a need in the art for achieving high throughput in a coded 
MIMO system without having to individually code different frequency bins of the 
spatial subchannels. 



SUMMARY 

[1009] Aspects of the invention provide techniques for processing a data 
transmission at the transmitter and receiver of a MIMO system such that high 
performance (i.e., high throughput) is achieved without the need to individually 
code/modulate for different frequency bins. In an aspect, a time-domain 
implementation is provided herein which uses frequency-domain singular value 
decomposition and "water-pouring" results to derive pulse-shaping and beam-steering 
solutions at the transmitter and receiver. The singular value decomposition is 
performed at the transmitter to determine the eigen-modes (i.e., the spatial subchannels) 
of the MIMO channel and to derive a first set of steering vectors that are used to 
"precondition" modulation symbols. The singular value decomposition is also 
performed at the receiver to derive a second set of steering vectors that are used to 
precondition the received signals such that orthogonal symbol streams are recovered at 
the receiver, which can simplify the receiver processing. Water-pouring analysis is used 
to more optimally allocate the total available transmit power for the MIMO system to 
the eigen-modes of the MIMO channel. The allocated transmit power may then 
determine the data rate and the coding and modulation scheme to be used for each 
eigen-mode. 

[1010] At the transmitter, data is initially coded in accordance with one or more 
coding schemes to provide coded data, which is then modulated in accordance with one 
or more modulation schemes to provide a number of modulation symbol streams (e.g., 
one stream for each eigen-mode). An estimated channel response matrix for the MIMO 
channel is determined (e.g., at the receiver and sent to the transmitter) and decomposed 
(e.g., in the frequency domain, using singular value decomposition) to obtain a first 
sequence of matrices of (right) eigen-vectors and a second sequence of matrices of 
singular values. Water-pouring analysis may be performed based on the matrices of 
singular values to derive a third sequence of matrices of values indicative of the transmit 
power allocated to the eigen-modes of the MIMO channel. A pulse-shaping matrix for 
the transmitter is then derived based on the first and third sequences of matrices. The 
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pulse-shaping matrix comprises the steering vectors that are used to precondition the 
modulation symbol streams to obtain a number of preconditioned signals, which are 
then transmitted over the MIMO channel to the receiver. 

[101 1) At the receiver, the estimated channel response matrix is also determined and 
decomposed to obtain a fourth sequence of matrices of (left) eigen-vectors, which are 
then used to derive a pulse-shaping matrix for the receiver. A number of signals is 
received at the receiver and preconditioned based on this pulse-shaping matrix to obtain 
a number of received symbol streams. Each received symbol stream may be equalized 
to obtain a corresponding recovered symbol stream, which is then demodulated and 
decoded to recover the transmitted data. 

[1012] Various aspects and embodiments of the invention are described in further 
detail below. The invention further provides methods, digital signal processors, 
transmitter and receiver units, and other apparatuses and elements that implement 
various aspects, embodiments, and features of the invention, as described in further 
detail below. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[1013] The features, nature, and advantages of the present invention will become 
more apparent from the detailed description set forth below when taken in conjunction 
with the drawings in which like reference characters identify correspondingly 
throughout and wherein: 

[1014] FIG. 1 is a block diagram of an embodiment of a transmitter system and a 
receiver system of a MIMO system; 

[1015] FIG. 2 is a block diagram of an embodiment of a transmitter unit capable of 

implementing various aspects and embodiments of the invention; 

[1016] FIG. 3 is a block diagram of an embodiment of a receiver unit capable of 

implementing various aspects and embodiments of the invention; 

[1017] FIGS. 4A and 4B are block diagrams of an equivalent channel model and an 

implementation, respectively, of a minimum mean square error linear equalizer 

(MMSE-LE); and 

[1018] FIGS. 5A and 5B are block diagrams of an equivalent channel model and an 
implementation, respectively, of a decision feedback equalizer (DFE). 
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DETAILED DESCRIPTION 
(1019) The techniques described herein for processing a data transmission at the 
transmitter and receiver may be used for various wireless communication systems. For 
clarity, various aspects and embodiments of the invention are described specifically for 
a multiple-input multiple-output (MIMO) communication system. 
(1020] A MIMO system employs multiple (Nt) transmit antennas and multiple (Nr) 
receive antennas for data transmission. A MIMO channel formed by the Nr transmit 
and Nr receive antennas may be decomposed into Nc independent channels, with Nc ^ 
min {Nt, Nr}. Each of the Nc independent channels is also referred to as a spatial 
subchannel (or a transmission channel) of the MIMO channel. The number of spatial 
subchannels is determined by the number of eigen-modes for the MIMO channel, which 
in turn is dependent on a channel response matrix that describes the response between 
the N T transmit and Nr receive antennas. 

[1021] FIG. 1 is a block diagram of an embodiment of a transmitter system 110 and 
a receiver system 150, which are capable of implementing various aspects and 
embodiments of the invention. 

[1022] At transmitter system 110, traffic data is provided from a data source 1 12 to 
a transmit (TX) data processor 114, which formats, codes, and interleaves the traffic 
data based on one or more coding schemes to provide coded data. The coded traffic 
data may then be multiplexed with pilot data using, e.g., time division multiplex (TDM) 
or code division multiplex (CDM) in all or a subset of the data streams to be 
transmitted. The pilot data is typically a known data pattern processed in a known 
manner, if at all. The multiplexed pilot and coded traffic data is then modulated (i.e., 
symbol mapped) based on one or more modulation schemes to provide modulation 
symbols, one modulation symbol stream for each spatial subchannel to be used for data 
transmission. The data rate, coding, interleaving, and modulation for each spatial 
subchannel may be determined by controls provided by a controller 130. 
[1023] The modulation symbols are then provided to a TX MIMO processor 120 
and further processed. In a specific embodiment, the processing by TX MIMO 
processor 120 includes (1) determining an estimated channel frequency response matrix 
for the MIMO channel, (2) decomposing the estimated channel frequency response 
matrix to determine the eigen-modes of the MIMO channel and to derive a set of 
"steering" vectors for the transmitter, one vector for the modulation symbol stream to be 
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transmitted on each spatial subchannel, (3) deriving a transmit spatio-temporal pulse- 
shaping matrix based on the steering vectors and a diagonal matrix indicative of the 
energies (i.e., the transmit powers) assigned to the eigen-modes, and (4) preconditioning 
(e.g., convolving) the modulation symbols with the pulse-shaping matrix to derive 
preconditioned modulation symbols. The processing by TX MIMO processor 120 is 
described in further detail below. Up to Nt streams of preconditioned modulation 
symbols are then provided to transmitters (TMTR) 122a through 122t. 
[1024] Each transmitter 122 converts the received preconditioned modulation 
symbol stream into one or more analog signals and further conditions (e.g., amplifies, 
filters, and quadrature modulates) the analog signals to generate a modulated signal 
suitable for transmission over the MEMO channel. The modulated signal from each 
transmitter 122 is then transmitted via a respective antenna 124 to the receiver system. 
[1025] At receiver system 150, the transmitted modulated signals are received by Nr 
antennas 152a through 152r, and the received signal from each antenna 152 is provided 
to a respective receiver (RCVR) 154. Each receiver 154 conditions (e.g., filters, 
amplifies, and downconverts) the received signal and digitizes the conditioned signal to 
provide a respective stream of samples. A RX MIMO processor 1 60 then receives and 
processes the sample streams to provide Nj streams of recovered modulation 
symbols. In an embodiment, the processing by RX MIMO processor 160 may include 

(1) determining an estimated channel frequency response matrix for the MIMO channel, 

(2) decomposing the estimated channel frequency response matrix to derive a set of 
steering vectors for the receiver, (3) deriving a receive spatio-temporal pulse-shaping 
matrix based on the steering vectors, (4) preconditioning (e.g., convolving) the samples 
with the pulse-shaping matrix to derive received modulation symbols, and (5) 
equalizing the received modulation symbols to derive the recovered modulation 
symbols. The processing by RX MIMO processor 160 is described in further detail 
below. 

[1026] A receive (RX) data processor 162 then demodulates, deinterleaves, and 
decodes the recovered modulation symbols to recover the transmitted traffic data. The 
processing by RX MIMO processor 160 and RX data processor 162 is complementary 
to that performed by TX MIMO processor 120 and TX data processor 1 14, respectively, 
at transmitter system 1 10. 
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(1027J RX MIMO processor 160 may further derive channel impulse responses for 
the MIMO channel, the signal-to-noise-and-interference ratios (SNRs) for the spatial 
subchannels, and so on, and provide these to a controller 170. RX data processor 162 
may also provide the status of each received packet or frame, one or more other 
performance metrics indicative of the decoded results, and possible other information. 
Controller 170 then derives channel state information (CSI), which may comprise all or 
some of the information received from RX MIMO processor 160 and RX data processor 
162. The CSI is processed by a TX data processor 178, modulated by a modulator 180, 
conditioned by transmitters 154a through 154r, and transmitted back to transmitter 
system 1 10. 

[1028] At transmitter system 110, the modulated signals from receiver system 150 
are received by antennas 124, conditioned by receivers 122, and demodulated by a 
demodulator 140 to recover the CSI transmitted by the receiver system. The CSI is then 
provided to controller 130 and used to generate various controls for TX data processor 
1 14 and TX MIMO processor 120. 

[1029] Controllers 130 and 170 direct the operation at the transmitter and receiver 
systems, respectively. Memories 132 and 172 provide storage for program codes and 
data used by controllers 130 and 170, respectively. 

[1030] In a MIMO system with limited total transmit power and operating over a 
frequency-selective channel (i.e., different gains at different frequencies), the channel 
capacity, C, is given by: 



where E T is the total available transmit power for the MIMO system; 

is an x ff R power spectral density matrix of an N R xl noise process 
vector, z(n) , at the receiver at frequency f k ; 
¥L(k) is an N K x N T channel frequency response matrix at frequency f k ; and 



C = max Slog | I+*nWH(*)*aWH W (*) I , 




Eq(l) 



subject to 



7 > 
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*„(£) is an N T xN T power spectral density matrix of an N T xl transmitted 
signal vector, x(w) , at frequency f k . 

[1031] The singular value decomposition (SVD) of the channel frequency response 
matrix H(k) at frequency f k may be expressed as: 



where U(/r) is an N R xN R unitary matrix (i.e., U // U=I, where I is the identity 
matrix with ones along the diagonal and zeros everywhere else); 
X(k) is an N R x N T diagonal matrix of singular values of H(fc) ; and 
V(£) is an N T x N T unitary matrix. 

The diagonal matrix \(k) contains non-negative real values along the diagonal (i.e., 
A.(£) = diag (^(k), ^(#), ... An T (k))) and zeros elsewhere. The A.(k) are referred to 
as the singular values of the matrix H(/c) . The singular value decomposition is a matrix 
operation known in the art and described in various references. One such reference is a 
book by Gilbert Strang entitled "Linear Algebra and Its Applications," Second Edition, 
Academic Press, 1980, which is incorporated herein by reference. 
[1032] It can be shown that in the uncorrelated white noise case (i.e., when 
N 

*„(£) = — -I , where No is the power spectral density of the noise at the receiver and 



l/T 0 is the bandwidth, in Hertz, of a frequency bin), the channel capacity is achieved 
when the power spectral density matrix *„(&) of the transmitted signal vector \(n) 
satisfies the condition: 



where E A (£) is an N T xN T diagonal matrix containing a set of energies (or transmit 

powers) assigned to the eigen-modes at frequency f k . The diagonal matrix E A (k) is 

the solution to the well-known "water-pouring" transmit energy distribution technique, 
which may be expressed as: 



Eq(2) 



*»(*) = ¥(*)!, (*)V ff (*) , 



Eq(3) 
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and 



Eq (4a) 



mjn(V-,/V r ) 



Eq (4b) 



where B is a constant derived from various system parameters. 

[1033] The water-pouring technique is analogous to pouring a fixed amount of 
water into a vessel with an irregular bottom, where each eigen-mode of each frequency 
bin corresponds to a point on the bottom of the vessel, and the elevation of the bottom at 
any given point corresponds to the inverse of the SNR associated with that eigen-mode. 
A low elevation thus corresponds to a high SNR and, conversely, a high elevation 
corresponds to a low SNR. The total available transmit power, Ej 9 is then "poured" into 
the vessel such that the lower points in the vessel (i.e., higher SNRs) are filled first, and 
the higher points (i.e., lower SNRs) are filled later. The constant B is indicative of the 
water surface level for the vessel after all of the total available transmit power has been 
poured, and may be estimated initially based on various system parameters. The 
transmit power distribution is dependent on the total available transmit power and the 
depth of the vessel over the bottom surface, and the points with elevations above the 
water surface level are not filled (i.e., eigen-modes with SNRs below a particular 
threshold are not used). 

[1034] The water-pouring technique is described by Robert G. Gallager, in 
"Information Theory and Reliable Communication," John Wiley and Sons, 1968, which 
is incorporated herein by reference. A specific algorithm for performing the basic 
water-pouring process for a MIMO-OFDM system is described in U.S. Patent 
Application Serial No. 09/978,337, entitled "Method and Apparatus for Determining 
Power Allocation in a MIMO Communication System," filed October 15, 2001, 
assigned to the assignee of the present application and incorporated herein by reference. 
[1035] The formulation of channel capacity shown in equations (1) through (4) 
suggests that an OFDM-based implementation may be used to achieve channel capacity 
by performing wateT-pouring in the frequency domain. With frequency-domain water- 
pouring, the total available transmit power is allocated to the N F frequency subchannels 
(or bins) on a bin-by-bin basis, with more power being allocated to bins achieving 
higher SNRs and less or no power being allocated to bins achieving lower SNRs. This 
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would then necessitate the use of separate coding and/or modulation schemes for each 
bin, which can complicate the coding and decoding at the transmitter and receiver, 
respectively. 

[1036] Aspects of the invention provide techniques for achieving high performance 
(i.e., channel capacity) via a time-domain implementation that uses frequency-domain 
singular value decomposition and water-pouring results to derive time-domain pulse- 
shaping and beam-steering solutions at the transmitter and receiver. 
[1037] The singular value decomposition is performed at the transmitter to 
determine the eigen-modes of the MEMO channel and to derive a first set of steering 
vectors that are used to precondition the modulation symbols. The singular value 
decomposition is also performed at the receiver to derive a second set of steering 
vectors that are used to precondition the received signals such that orthogonal symbol 
streams are recovered at the receiver, which can simplify the receiver processing. 
Water-pouring analysis is used to more optimally allocate the total available transmit 
power for the MIMO system to the eigen-modes such that high performance is 
achieved. The allocated transmit power may then determine the data rate and the 
coding and modulation scheme to be used for each eigen-mode. 

[1038] The techniques described herein provide several potential advantages. First, 
with time-domain eigen-mode decomposition, the maximum number of data streams 
with different SNRs, and thus different coding/modulation requirements, is given by 
min(A^ r ,^). It is also possible to make the received SNRs for the data streams 

essentially the same, thereby further simplifying the coding/modulation. The inventive 
techniques can thus greatly simplify the coding/modulation for a data transmission by 
avoiding the per-bin bit allocation required to approach channel capacity in OFDM 
systems that utilize frequency-domain water-pouring. 

[1039] Second, the orthogonalization process at the receiver results in decoupled 
(i.e., orthogonal) received symbol streams. This then greatly reduces the complexity of 
the time-domain equalization required for the decoupled symbol streams. In this case, 
the equalization may be achieved by simple linear space-time processing followed by 
parallel time-domain equalization of the independent symbol streams. In contrast, other 
wide-band time-domain techniques typically require more complicated space-time 
equalization to recover the symbol streams. 
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[10401 Third, the time-domain signaling techniques of the invention can more easily 
integrate the channel/pilot structures of various CDMA standards, which are also based 
on time-domain signaling. Implementation of channel/pilot structures may be more 
complicated in OFDM systems that perform frequency-domain signaling. 
[1041] FIG. 2 is a block diagram of an embodiment of a transmitter unit 200 
capable of implementing various aspects and embodiments of the invention. 
Transmitter unit 200 is an embodiment of the transmitter portion of transmitter system 
1 1 0 in FIG. 1 . Transmitter unit 200 includes (1) a TX data processor 1 14a that receives 
and processes traffic and pilot data to provide N T modulation symbol streams and (2) a 
TX MIMO processor 120a that preconditions the modulation symbol streams to provide 
^preconditioned modulation symbol streams. TX data processor 1 14a and TX MIMO 
processor 120a are one embodiment of TX data processor 1 14 and TX MIMO processor 
120, respectively, in FIG. 1. 

[1042] In the specific embodiment shown in FIG. 2, TX data processor 114a 
includes an encoder 212, a channel interle aver 214, and a symbol mapping element 216. 
Encoder 212 receives and codes the traffic data (i.e., the information bits, 6 f ) in 

accordance with one or more coding schemes to provide coded bits. The coding 
increases the reliability of the data transmission. In an embodiment, a separate coding 
scheme may be used for the information bits for each spatial subchannel. In alternative 
embodiments, a separate coding scheme may be used for each subset of spatial 
subchannels, or a common coding scheme may be used for all spatial subchannels. The 
coding scherne(s) to be used is determined by the controls from controller 130, which 
may be determined based on CSI received from the receiver system. Each selected 
coding scheme may include any combination of cyclic redundancy check (CRC), 
convolutional coding, Turbo coding, block coding, and other coding, or no coding at all. 
[1043] Channel interleaver 214 interleaves the coded bits based on one or more 
interleaving schemes (e.g., one interleaving scheme for each selected coding scheme). 
The interleaving provides time diversity for the coded bits, permits the data to be 
transmitted based on an average SNR for each spatial subchannel used for the data 
transmission, combats fading, and further removes correlation between coded bits used 
to form each modulation symbol. 

[1044] Symbol mapping element 216 then receives and multiplexes pilot data with 
the interleaved data and further maps the multiplexed data in accordance with one or 
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more modulation schemes to provide modulation symbols. A separate modulation 
scheme may be used for each spatial subchannel or for each subset of spatial 
subchannels. Alternatively, a common modulation scheme may be used for all spatial 
subchannels. The symbol mapping for each spatial subchannel may be achieved by 
grouping sets of bits to form non-binary symbols and mapping each non-binary symbol 
to a point in a signal constellation corresponding to the modulation scheme (e.g., QPSK, 
M-PSK, M-QAM, or some other scheme) selected for that spatial subchannel. Each 
mapped signal point corresponds to a modulation symbol. Symbol mapping element 
216 provides a vector of modulation symbols for each symbol period, with the number 
of modulation symbols in each vector corresponding to the number of spatial 
subchannels selected for use for that symbol period. Symbol mapping element 216 thus 
provides up to Nt modulation symbol streams (i.e., a sequence of symbol vectors with 
each vector including up to Nt modulation symbols), which are also referred to herein as 
the transmitted symbol vector, s(n) . 

[1045] The response of the MIMO channel to be used for data transmission is 
estimated and used to precondition the transmitted symbol vector prior to transmission 
to the receiver system. In a frequency division duplex (FDD) system, the downlink and 
uplink are allocated different frequency bands, and the responses for the downlink and 
uplink may not be correlated to a sufficient degree. For the FDD system, the channel 
response may be estimated at the receiver and sent back to the transmitter. In a time 
division duplex (TDD) system, the downlink and uplink share the same frequency band 
in a time division multiplexed manner, and a high degree of correlation may exist 
between the downlink and uplink responses. For the TDD system, the transmitter 
system can estimate the uplink channel response (e.g., based on the pilot transmitted by 
the receiver system on the uplink) and derive the downlink channel response by 
accounting for differences between the transmit and receive antenna array manifolds. 
[1046] In an embodiment, channel response estimates are provided to TX MIMO 
processor 120a as a sequence of N R x N T matrices of time-domain samples, &_(ri) . 
The (/, j) -th element of the estimated channel impulse response matrix 0£_(n) , for 
l<,i^N R and 1 £ j <, N T , is a sequence of samples representing the sampled impulse 
response of the propagation path from the y-th transmit antenna to the i-th receive 
antenna. 
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[1047] Within TX MIMO processor 120a, a fast Fourier transformer 222 receives 
the estimated channel impulse response matrix, 2L(n) , (e.g., from the receiver system) 
and derives the corresponding estimated channel frequency response matrix, H(fc) , by 

performing a fast Fourier transform (FFT) on jjT(n) (i.e., H(Jfc) = FFT[#(/i)]). This 
may be achieved by performing an Afc-point FFT on a sequence of N F samples for each 
element of 2[{n) to derive a sequence of Nf coefficients for the corresponding element 
of H(fr) . The Nr-N t elements of H{k) are thus Nr-Nt sequences representing the 
frequency responses of the propagation paths between the N T transmit antennas and Nr 
receive antennas. Each element of H(/r) is the FFT of the corresponding element of 

[1048] A block 224 then computes the singular value decomposition of the 
estimated channel frequency response matrix, H(£), for each value of k> where 
0 ^ k £ (N F -1) and Nf is the length of the FFT (i.e., N F corresponds to the number of 
frequency bins). The singular value decomposition may be expressed as shown in 
equation (2), which is: 

H(k) = V(k)\(k)V H (k) . 

The result of the singular value decomposition is three sequences of N? matrices, 
U(ifc), \(k), and V*(*), for 0£k£(N F -l). For each value of k, V(k) is the 
N R xN R unitary matrix of left ei gen-vectors of H(fc), V(fc) is the N T xN T unitary 
matrix of right eigen-vectors of H(£), and \(k) is the N R xN T diagonal matrix of 
singular values of H(k) . 

[1049] The singular value decomposition is used to decompose the MIMO channel 
into its eigen-modes, at the frequency f k associated with frequency bin k , for each 

value of 0 ^ k £ (N F -\) . The rank r(k) of H(/r) corresponds to the number of 
eigen-modes for the MIMO channel at frequency f k , which corresponds to the number 
of independent channels (i.e., the number of spatial subchannels) available in frequency 
bin k. As described in fiirther detail below, the columns of V(fc) are the steering 
vectors associated with frequency f k to be used at the transmitter for the elements of 
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the transmitted symbol vector, s(n) . Correspondingly, the columns of V(k) are the 
steering vectors associated with frequency f k to be used at the receiver for the elements 
of the received signal vector, r(/i). The matrices U(&) and V(A), for 
OZk <(N F -l), are used to orthogonalize the symbol streams transmitted on the 
eigen-rnodes at each frequency f k . When these matrices are used collectively to 

preprocess the transmitted and received symbol streams, either in the frequency domain 
or the time domain, as described in detail below, the result is the overall 
orthogonalization of the received symbol stream. This then allows for separate 
coding/modulation per eigen-rnode (as opposed to per bin) and further simplification of 
the equalization of the received symbol streams at the receiver, as described below. 
[1050] The elements along the diagonal of X(k) are A l7 (£) for l£i < r(k) , where 

r{k) is the rank of H(Jt). The columns of U(fc) and V(*), u,(fc) and v,.(A;), 
respectively, are solutions to the eigen equation, which may be expressed as: 

fi(*)v,(*)-^a(*) • Eq(5) 

[1051] The A(£), and V(#) matrices may be provided in two forms - a 

"sorted" form and a "random-ordered" form. In the sorted form, the diagonal elements 
of \(k) are sorted in decreasing order so that A^ y (k)> A^ 2 (k)>K >A rr (k) 9 and their 
eigen-vectors are arranged in corresponding order in V(k) and V(*) . The sorted form 
is indicated herein by the subscript s y i.e., V s (k) 9 h s (k), and V s (k). In the random- 
ordered form, the ordering of the singular values and eigen-vectors is random and 
independent of frequency. The random form is indicated herein by the subscript r. The 
particular form selected for use, sorted or random-ordered, determines the eigen-modes 
to be used for the data transmission and the coding and modulation scheme to be used 
for each selected eigen-mode. 

[1052] A water-pouring analysis block 226 then receives the set of singular values 
for each frequency bin, which is contained in the sequence of matrices, X(fc), and CSI 
that includes the received SNR corresponding to each singular value. The received 
SNR is the SNR achieved at the receiver for the recovered modulation symbols, as 
described below. The matrices \(k) are used in conjunction with the received SNRs to 
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derive the sequence of diagonal matrices, E x {k) 9 which are the solution to the water- 
pouring equations (4a) and (4b). As noted above, the diagonal matrices E x (k) contain 
the set of energies or transmit powers assigned to the ei gen-modes at each of the N F 
frequency bins. The water-pouring analysis used to derive the diagonal matrices, 
E x (k), may be performed as described in the aforementioned U.S. Patent Application 
Serial No. [Attorney Docket No. 010467]. 

[1053] A scaler/ffFT 228 receives the unitary matrices, V(fc), and the diagonal 
matrices, (A), for all Nf frequency bins, and derives a spatio-temporal pulse-shaping 
matrix, P^rt) , for the transmitter based on the received matrices. Initially, the square 
root of the diagonal matrices, E x (k) y is computed to derive a sequence of diagonal 
matrices, ^JE X (k) , whose elements are the square roots of the elements of E x (k). The 
elements of the diagonal matrices, E^(A), are representative of the transmit power 
allocated to the eigen-modes. The square root then transforms the power allocation to 
the equivalent signal scaling. The product of the square-root diagonal matrices, 

yjE x (k) , and the unitary matrices, V(k) , which are the sequence of matrices of right 
eigen-vectors of H(£) , is then computed. This product, V(k)yjE x (k) , defines the 
optimal spatio-spectral shaping to be applied to the transmitted symbol vector, s(n) . 
[1054] An inverse FFT of the product V(k)^]E x (k) is then computed to derive the 
spatio-temporal pulse-shaping matrix, P^A), for the transmitter, which may be 
expressed as: 

P tf (X) = IFFT[y(*)VE^ft)] - Eq(6) 

The pulse-shaping matrix, P te (A) , is an N T xN T matrix. Each element of P tt (A) is a 
sequence of values. Each column of P^A.) is a steering vector for a corresponding 
element of s(«) . 

[1055] A convolver 230 receives and preconditions (e.g., convolves) the transmitted 
symbol vector, s(n) 9 with the pulse-shaping matrix, P tt (A), to derive the transmitted 
signal vector, x(/a) . The convolution of s(w) with P a (A) may be expressed as: 
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Eq(7) 



The matrix convolution shown in equation (7) may be performed as follows. To derive 
the Mh element of the vector x(n) for time n, the inner product of the z-th row of 

the matrix P tt (X) with the vector s{n - A,) is formed for a number of delay indices (e.g., 
0< X< (N F -1)), and the results are accumulated to derive the element *,.(«). The 
signal transmitted on each transmit antenna (i.e., each element of x(n) or *, («)) is thus 
formed as a weighted combination of the Nr modulation symbol streams, with the 
weighting determined by the appropriate column of the matrix P tt (k) ■ The process is 
repeated such that each element of the vector x(n) is derived from a respective column 
of the matrix P tt (X) and the vector s(/i) . 

[1056] Each element of the transmitted signal vector, x(n) , corresponds to a 
sequence of preconditioned symbols to be transmitted over a respective transmit 
antenna. The N T preconditioned symbol sequences (i.e., a sequence of preconditioned 
symbol vectors with each vector including up to Nt preconditioned symbols) correspond 
to Nr transmitted signals, and are also referred to herein as the transmitted signal vector, 
x(n) . The N T transmitted signals are provided to transmitters 122a through 122t and 
processed to derive Nt modulated signals, which are then transmitted from antennas 
124a through 124t, respectively. 

[1057] The embodiment shown in FIG. 2 performs time-domain beam-steering of 
the transmitted symbol vector, s(n) . The beam-steering may also be performed in the 
frequency domain. In this case, the vector s(/i) may be transformed via an FFT to 
derive a frequency-domain vector S(k) . The vector S(k) is then multiplied with the 
matrix y(k)yjE^Xk) to derive a frequency-domain vector X(k) , as follows: 

The transmitted signal vector, x(n) , may then be derived by performing an IFFT on the 
vector X(ifc) (i.e., x(n) = IFFT[X(*)]). 

[1058] FIG. 3 is a block diagram of an embodiment of a receiver unit 300 capable of 
implementing various aspects and embodiments of the invention. Receiver unit 300 is 
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an embodiment of the receiver portion of receiver system 150 in FIG. 1- Receiver unit 
300 includes (1) a RX MIMO processor 1 60a that processes Nr received sample streams 
to provide N T recovered symbol streams and (2) a RX data processor 162a that 
demodulates, deinterleaves, and decodes the recovered symbols to provide decoded bits. 
RX MIMO processor 160a and RX data processor 162a are one embodiment of RX 
MIMO processor 160 and RX data processor 162, respectively, in FIG. 1. 
[1059] Referring back to FIG. 1, the transmitted signals from iv> transmit antennas 
are received by each of antennas 152a through 152r, and the received signal from 
each antenna is routed to a respective receiver 154 (which is also referred to as a front- 
end processor). Each receiver 1 54 conditions (e.g., filters and amplifies) a respective 
received signal, downconverts the conditioned signal to an intermediate frequency or 
baseband, and digitizes the downconverted signal to provide ADC samples. Each 
receiver 154 may further data demodulate the ADC samples with a recovered pilot to 
generate a respective stream of received samples. Receivers 154a through 154r thus 
collectively provide Nr received sample streams (i.e., a sequence of vectors with each 
vector including up to Nr samples), which are also referred to as the received signal 
vector, r(n) . The received signal vector, r(n) , is then provided to RX MIMO 
processor 160a. 

[1060] Within RX MEMO processor 160a, a channel estimator 312 receives the 
vector r(n) and derives an estimated channel impulse response matrix, ti(n) , which 
may be sent back to the transmitter system and used in the transmit processing. An FFT 
314 then performs an FFT on the estimated channel impulse response matrix:, 2£(n) , to 
derive the estimated channel frequency response matrix, U(k) . 

[1061] A block 316 then computes the singular value decomposition of H(£) for 
each value of k to obtain the matrix of left eigen-vectors, V(k) , for the corresponding 
frequency bin k. Each column of V(k) is a steering vector for a corresponding 
element of r(«), and is used to orthogonalize the received symbol streams at the 
receiver system. An IFFT 318 then performs the inverse FFT of V(k) to derive a 
spatio-temporal pulse-shaping matrix, V(X) , for the receiver system. 
[1062] A convolver 320 then derives the received symbol vector, r(n) , which is an 
estimate of the transmitted symbol vector, s(n) , by performing a convolution of the 
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received signal vector, r(n) , with the conjugate transpose of the spatio-temporal pulse- 
shaping matrix, V H (\). This convolution may be expressed as: 



[1063] The pulse-shaping at the receiver may also be performed in the frequency 
domain, similar to that described above for the transmitter. In this case, the received 
signal vector, r(«), may be transformed via an FFT to derive a frequency-domain 
vector R(A) . The vector R(k) is then pre-multiplied with the conjugate transpose 
matrix V H {k) to derive a frequency-domain vector R(k) . The result of this matrix 
multiplication, R(fc), can then be transformed via an inverse FFT to derive the time- 
domain received symbol vector, r(n) . The convolution of the vector r(/i) with the 
matrix V H (A.) can thus be represented in the discrete frequency domain as: 



where \(k) = X(k)^E x {k) is a matrix of weighted singular values of U(k) , with the 



S(k) is the FFT of s(n) , the transmitted symbol vector; 
R(A) is the FFT of r(n) , the received signal vector; 
R(*) is the FFT of r(n) , the received symbol vector; 
Z(k) is the FFT of z(n) , the vector of received noise samples; and 
Z{k) is the FFT of the received noise process as transformed by the unitary 
matrix U w (k). 

[1064] From equation (9), the received symbol vector, r(/i), may be characterized 
as a convolution in the time domain, as follows: 



Eq(8) 



to = U"WEW = ^)S(/t) + Z(Ar) , 



Eq(9) 



weights being the square root of the water-pouring solution, -^E^ (k) ; 



f(") = £A(X)s(/>-X) + z(«) , 



Eq(10) 



where A (X) is the inverse FFT of X(k) = \(k)jEjk) ; and 
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i(n) is the received noise as transformed by the receiver spatio-temporal pulse- 
shaping matrix, V" (A). 
The matrix A(X) is a diagonal matrix of eigen-pulses, with each such eigen-pulse being 
derived as the EFFT of the corresponding set of singular values in \(k) for 
0£*£(tf r -l). 

[1065] The two forms for ordering the singular values, sorted and random-ordered, 
result in two different types of eigen-pulses. For the sorted form, the resulting eigen- 
pulse matrix, A s (1 ) , is a diagonal matrix of pulses that are sorted in descending order of 
energy content. The pulse corresponding to the first diagonal element of the eigen-pulse 
matrix, { A ,(1 )} n , has the most energy, and the pulses corresponding to elements further 

down the diagonal have successively less energy. Furthermore, when the SNR is low 
enough that water-pouring results in some of the frequency bins having no energy, the 
energy is taken out of the smallest eigen-pulses first. Thus, at low SNRs, one or more 
of the eigen-pulses may have no energy. This has the advantage that at low SNRs, the 
coding and modulation are simplified through the reduction in the number of orthogonal 
subchannels. However, in order to approach channel capacity, it is necessary to code 
and modulate separately for each eigen-pulse. 

[1066] The random-ordered form of the singular values in the frequency domain 
may be used to further simplify coding and modulation (i.e., to avoid the complexity of 
separate coding and modulation for each element of the eigen-pulse matrix). In the 
random-ordered form, for each frequency bin, the ordering of the singular values is 
random instead of being based on their size. This random ordering can result in 
approximately equal energy in all of the eigen-pulses. When the SNR is low enough to 
result in frequency bins with no energy, these bins are spread approximately evenly 
among the eigen-modes so that the number of eigen-pulses with non-zero energy is the 
same independent of SNR. At high SNRs, the random-order form has the advantage 
that all of the eigen-pulses have approximately equal energy, in which case separate 
coding and modulation for different eigen-modes is not required. 
[1067] If the response of the MIMO channel is frequency selective (i.e., different 
values for H(A) for different values of A), the eigen-pulses in the matrix A (A) are time- 
dispersive. In this case, the resulting received symbol sequences, r(n) , have inter- 
symbol interference (IS I) that will in general require equalization to provide high 
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performance. Furthermore, because the singular values in X(k) are real, the elements of 
\(k) = X(k)ylE A (k) are also real, and the eigen-pulses in the matrix A(l) exhibit 

aliased conjugate symmetry properties. If steps are taken to avoid this time-domain 
aliasing (e.g., by using an FFT length, N Fy that is sufficiently greater than the number of 

non-zero samples in the estimated channel impulse response matrix, 2£(n) ) then the 
eigen-pulse matrix is conjugate symmetric in the delay variable, 1 , i.e., A (A.) = A*(-A) . 
[1068] An equalizer 322 receives the received symbol vector, r(w) , and performs 
space-time equalization to derive a recovered symbol vector, s(/i) , which is an estimate 
of the transmitted symbol vector, s(n) . The equalization is described in further detail 
below. The recovered symbol vector, s(n) , is then provided to RX data processor 162a. 
[1069] Within RX data processor 162a, a symbol unmapping element 332 
demodulates each recovered symbol in s(«) in accordance with a demodulation scheme 
(e.g., M-PSK, M-QAM) that is complementary to the modulation scheme used for that 
symbol at the transmitter system. The demodulated data from symbol unmapping 
element 332 is then de-interleaved by a deinterleaver 334, and the deinterleaved data is 

further decoded by a decoder 336 to obtain the decoded bits, b t , which are estimates of 

the transmitted information bits, b t . The deinterleaving and decoding are performed in 

a manner complementary to the interleaving and encoding, respectively, performed at 
the transmitter system. For example, a Turbo decoder or a Viterbi decoder may be used 
for decoder 336 if Turbo or convolutional coding, respectively, is performed at the 
transmitter system. 



Minimum Mean Square Error (MMSE) Equalization 

(1070] As shown in equation (10), an equivalent channel for the received symbol 
vector, r(n) , has an impulse response (i.e., a unit sample response) of A(X) , which is 
the diagonal matrix of eigen-pulses, and a corresponding frequency response of X(/) . 
A matched filter receiver for r(n) would then include a filter matched to the impulse 
response of A(X) . Such a matched filter would have an impulse response of A"(-l) 
and a frequency response of X'(/) , which may be expressed as: 
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XV^^A"*!)^ • Eq(ll) 

The end-to-end frequency response of the equivalent channel for f(«) and its matched 
filter may be given as ^(/) = A(/)A'(/) • 

[1071] The end-to-end frequency response of \p(f) may be spectrally factorized 

into a hypothetical filter and its matched filter. This hypothetical filter would have a 
causal impulse response of £(X), where £(X) = 0 for A<0 , and a frequency response 
of y{f) . The end-to-end frequency response of the hypothetical filter and its matched 
filter is (by definition) equal to the end-to-end frequency response of the equivalent 
channel and its matched filter, i.e., y(f)y H (f) = . 

[1072] For the following analysis, an equivalent channel model may be defined to 
have spectrally white noise. This may be achieved by applying a noise-whitening filter 

having a frequency response matrix, (V'C/)) = {y(f)y H {/)) y(f)> which is the 
Moore-Penrose inverse of y H (f), to the output of the receiver matched filter. The 
overall frequency response of the channel (with frequency , response of A(/)), the 
matched filter (with frequency response of X'(/) ), and the noise- whitening filter (with 

frequency response of (7" (/)) ) may then be expressed as: 

A(/)X'(/)(f(/)f = *(/)(/ (/)) + = y(f) . Eq (12) 

The impulse response £(X) corresponding to the frequency response y(f) is a diagonal 
matrix. 

[1073] FIG- 4A is a diagram of a minimum mean square error linear equalizer 
(MMSE-LE) 414 derived based on an equivalent channel model. The received symbol 
vector, r(n), is filtered by a (hypothetical) whitened matched filter 412 to provide a 
filtered symbol vector, ?(/i). Whitened matched filter 412 performs the dual function 
of matched filtering for r(n) and noise whitening, and has a response of 
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X'(/) {j" (/))*• The filtered symbol vector, ?(«), is the output of the equivalent 
channel model and may be expressed as: 



l(n) = J] £(X) s(n - X) + z(/i) = £ s(«) + z(n) , 



Eq(13) 



where £ is an N R x(L + \)N T block-structured matrix that represents the sequence of 

matrices, , for the sampled channel-whitened eigen-pulses and can be represented 

as: 

£ = [£(0) E(l) A £(£)] , 
and s(«) is a sequence of £+1 vectors of modulation symbols and can be represented as: 



S(») 
§(«-!) 

M 
s(n-L) 



Each vector of s(n) comprises up to Nt symbols and each symbol in the vector is 
associated with one of the eigen-pulses in the matrix V . The blocks of £ (i.e., 
E(0), C(l), A , £(£)) are all diagonal. 

[1074) When the receiver input noise is white and has a power spectral density of 
Af 0 I, the noise vector i(n) has an autocorrelation function, q>„ (Ar) , which can be 

expressed as: 



%,(m) = E[i(n-m)i H (n)] = N 0 ^ vv (m) , 



Eq(14) 



where 
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Since the sequence of matrices, V(&) , of right eigen-vectors of H(fc) are all unitary, 
then V"(fc)V(fc) = I for each value of k. As a result, <p w (/n) , which is the inverse FFT 
of the sequence V W (A:)V(/:) , is given by: 



where 5{m) is the unit sample sequence, which may be expressed as: 



Eq(15) 



S(m) = 



(1, w = 0 
1 0 , otherwise 



The noise vector, z(«), after the whitened matched filter has an autocorrelation 



function, <P~("0 > which may then be expressed as: 



<P^(m) = 9M = N 0 lS(m) 



Eq(16) 



[1075] An MMSE-LE computes an initial estimate, , of the transmitted symbol 
vector, s(n) , at time n by performing a matrix convolution of the sequence of filtered 
symbol vectors, r(>i), with a sequence of 2K + 1 , M T xN R weight matrices, M(A), as 
follows: 



?(«) = £m(A.) ? (* = MEM , 



Eq(17) 



where M=[M(~^) A M(0) A M(tf)]; 

ATis a parameter that determines the delay-extent of the equalizer; and 



M 
M 



|1076] The sequence of weight matrices, M(X) , is selected to minimize the mean- 
square error, which can be expressed as: 
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s = E{e M (n)e(n)} , Eq(18) 
where the error e(n) can be expressed as: 

e(") = s(«)-s(*) . Eq(19) 

[1077] The MMSE solution can then be stated as the sequence of weight matrices, 
M(l ) , that satisfy the following linear constraints: 



XM(X)2_a-rfl)== 



0, -K<m<-L 
E"(-ro), -L^m^O 
0, 0<m£K 



E"(-ro), -L<S/w<S0 , Eq(20) 



where V~(rn) is a sequence of N R x N R space-time correlation matrices. The matrices 
(p^irn) can be expressed as: 



? ?? ("0 = E[L{n - rn)r_ H («)] = < 



£E(X)r(X+w)+2 (/w) , -L<>m<L 

, Eq(21) 

q> (/w) , otherwise 



where <P—(tf0 is given by equations (14) through (1 6). 

[1078] For spatially and temporally uncorrelated noise, <P~(w) = N Q \S(m) . In this 
case, all of the off-diagonal terms in 9-^-™) and r r (-/w) are zero, and all of the 
off-diagonal terms in M(A) are also zero, which then results in a set of decoupled 
equations for the equalizer coefficients. Thus, the linear constraints in equation (20) 
maybe simplified as follows: 

fc«irW(j ff (*.-*)\=n(-*) . forl<**<r; -L<mZ0 Eq (22) 

where r is the rank of the matrix H(/r) . 

[1079] Equation (20) can further be represented as: 

M<p =T\ or M = r // q>" 1 , Eq (23) 
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where <p_ is block-Toeplitz with block j,k given by <p~(j-k) and 



r = 



XL) 

£(i-D 

M 

£(0) 



—KN R xN T 



where 0 mx „ is an mxn matrix of zeros. 

[1080] The frequency response matrix, m(f) , corresponding to the time-domain 
weight matrices, M(X), —L<> \<L, for the MMSE-LE may be derived by taking the 
matrix Fourier transform of M (A) , as follows: 



Since M(A) is diagonal, the frequency response matrix, m(/) , is also diagonal. 
[1081] As shown in FIG. 4A, the filtered symbol vector, ?(«), is provided to an 
MMSE-LE 414 and equalized based on the frequency response matrix, m(/) , to derive 
the symbol vector, s(h), which is the estimate of the transmitted symbol vector, s(w) . 
Because of the pulse-shaping performed at both the transmitter and receiver systems, 
the received symbol sequences in r(n) are orthogonal and the weight matrices, M(X) , 
for the MMSE-LE are diagonal matrices. Thus, each of the Nr received symbol 
sequences in r(n) may be equalized independently by the MMSE-LE, which can 
greatly simplify the receiver processing. 

[1082] To determine the SNR associated with the symbol estimates, ?(w), an 
unbiased minimum mean square error estimate is first derived. For the initial symbol 
estimate, ?(/?) , derived above, 



L 



m(/)=£Mw e - 



Eq(24) 



^[l(«)|i(«)]=M£[r(«)|s(«)] 

= [M(-K)rs(/i + K) +A + M(0)rs(/j) +A + M(/0£s(« - * )] 



Eq(25) 



WO 03/050968 PCT/US02/38671 
020068 

26 

where the expectation is taken over the noise. If it is assumed that the modulation 
symbols are uncorrelated in time and the expectation is taken over all inter-symbol 
interference in the above (all transmitted signal components not transmitted at time n\ 
then the expectation can be expressed as: 

£[!(") I §(")] =M£[r(*)|s(*)] 

= [M(0)E(0) + M(-1)E(1)+ - +M(-L)E(£)]s(n) 

. Eq (26) 

= MTs(n) 
= Gs(«) 

where 

G = M£ = r /f 9^r . 

When the noise is spatially and temporally uncorrelated, M(A), for -K <>X£K , are 
diagonal, so G is N T x N T diagonal. 

[1083J After averaging over the interference from other spatial subchannels, the 
average value of the signal from the i-th transmit antenna at time n can be expressed as: 

£[?(») I *,(«)] , Eq(27) 

where g u is the i-th diagonal element of G (g u is a scalar), and s { (n) is the i-th 
element of the initial symbol estimate, ?(/i) . 
[1084] By defining 

D^ 1 =diag(l/ g|1 , l/g n9 A , \lg NrNj ) , Eq(28) 

the unbiased symbol estimate, s(n) , of the transmitted symbol vector, s(n) , at time n 
may then be expressed as: 

sW = Dci(«) = DcMr(«) . Eq(29) 



[1085] The error covariance matrix associated with the unbiased symbol estimate, 
s(w) , can be expressed as: 
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? Cf = W = £{[s(n)-D- I Mr(n)] [s H («)- fW^D"'] } 

. Eq (30) 

^-DcG-GD-'+D-'GD^ 

For the spatially and temporally uncorrelated noise case, =G"\ so in this case 
W = G 1 -I 

[1086] The SNR associated with the unbiased estimate, $.(n) , of the symbol 
transmitted on the i-th transmit antenna can finally be expressed as: 

SNR. = — = . Eq (31) 

[1087] In FIG. 4A, whitened matched filter 412 in the equivalent channel model is 
provided to simplify the derivation of the MMSE-LE. In a practical implementation, 
the response of the whitened matched filter is (automatically) incorporated within the 
response of the MMSE-LE when the MMSE-LE is adapted to minimize the mean 
square error. 

[1088] FIG. 4B is a block diagram of an embodiment of an MMSE-LE 322a, which 
is an embodiment of equalizer 322 in FIG. 3. Initially, the matrices H and q>_ may 

first be estimated based on the received pilot and/or data transmissions. The weight 
matrices M are then computed according to equation (23). 

[1089] Within MMSE-LE 322a, the received symbol vector, r(n) , from RX MIMO 
processor 160 are pre-multiplied by a multiplier 422 with the weight matrices M to 
form an initial estimate, s(/i) , of the transmitted symbol vector, s(w) , as shown above 
in equation ( 1 7). The initial estimate, l(n) , is further pre-multiplied by a multiplier 424 
with the diagonal matrix D^ 1 to form an unbiased estimate, s(w), of the transmitted 
symbol vector, s(n) , as shown above in equation (29). The unbiased estimate, s(/i) , 
comprises the recovered symbol vector that is provided by the MMSE-LE to RX data 
processor 162. 

[1090] The recovered symbol vector, s(n) 9 is also provided to a CSI processor 428, 
which derives CSI for the MIMO channel. For example, CSI processor 428 may 
estimate the SNR of the i-th recovered symbol sequence according to equation (31). 
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The SNRs for the recovered symbol sequences comprise a part of the CSI that is 
reported back to the transmitter unit. 

[1091] The recovered symbol vector, s(h), are further provided to an adaptive 
processor 426, which then derives the weight matrices M and the diagonal matrix 
based on equation (23) and (28), respectively. 



Decision Feedback Equalization 

[1092] A decision feedback equalizer (DFE) used in conjunction with the wideband 
eigen-mode transmission forms an initial estimate, s(w), of the transmitted symbol 
vector, s(/i) , at time n, which can be expressed as: 



100= ZM/Wit/i-^ + ^M.wk/i-X) t 



Eq(32) 



where f$*i) is the vector of filtered modulation symbols given by equation (13); 

k(/i) is a vector of remodulated symbols (i.e., symbols that have been 

demodulated and then modulated again); 
M f (X), -AT, <X^0, is a sequence of (tf, +1) - N T x N R feed-forward 

coefficient matrices; and 
M b (k), \<X£K 29 is a sequence of K 2 -N T *N R feed-back coefficient 
matrices. 

Equation (32) can also be expressed as: 



s(«) = M / r(/i) + M A i(n) , 

where M f = M f (-K x +1)L M ,(<))] ; 

M 6 =[M,(1)M,(2)L M b (K 2 )] ; 



{(n) 



Eq(33) 









i(n-2) 
M 


; and £fa) = 


M 
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[1093] If the MMSE criterion is used to determine the feed-forward and feed-back 
coefficient matrices, then the solutions for M and M _ that minimize the mean square 

error, 

can be used, where the error e(/i) is expressed as: 



e(rt) = - s(w) . 



Eq (34) 



[1094] The MMSE solution for the feed-forward filter, M f (X), for - K x £ X£ 0 , is 
determined by the following linear constraints: 

£ M^OrXECOr^/^l-^ + ^O^J^r^C-m) , Eq(35) 

I—*, L«=o J 



and can also be expressed as: 



M.=r"9_" . 

/ — =rr 



Eq (36) 



where 



r = 



E(L) 
E(L-l) 
M 

E(0) 



and 9_ is a (K l +l)N R x(K s +l)N K matrix made up of N R xN,, blocks. The (i,»-th 

= rr 

block in <p_ is given by: 



Ai-ITI 



Eq(37) 



X*0 



[1 095] The MMSE solution for the feed-back filter may be expressed as: 
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M 4 (X) = - Jm, (/)£(*.-/) , {or\<X<.K 2 



i—K, 



Eq(38) 



where f = [ff 2 L f L Q {Kl+mKl _ l)Nl ] ,and f x = 



£(L) 

M 
C(X) 



[1096] Since the matrices £(X), forO^X^Z, are diagonal, then from equation 
(36), the feed-forward filter coefficient matrices, My (A), for -AT, £A^0, are also 
diagonal. It then follows that the feed-back filter coefficient matrices, 
M b {X}, for 1 <, X< K 29 are also diagonal. 

[1097] The feed- forward filter and feed-back filter have frequency response 
matrices m f (f) and m 6 (/) , respectively, which are given by: 



m f (f) = lM ; (X)e- w ,and 
m 6 (/) = £mi(^" w • 



Eq(39) 



[1098] FIG. 5A is a diagram of a decision feedback equalizer derived based on an 
equivalent channel model. The received symbol vector, r(/i), is filtered by a 
(hypothetical) whitened matched filter 512 to provide the filtered symbol vector, f(/t). 
The vector r_(n) is further filtered by a feed-forward filter 514 having the frequency 
response of m f {f}. The output from feed-forward filter 514 is summed with the 
output from a feed-back filter 518 by a summer 516 to derive the symbol vector \(ri). 
This vector s(w) is also provided to a symbol decision element 520 to derive the 
remodulated symbol vector, ^(n) , which represents the detected symbols for the symbol 
estimate s(«). The remodulated symbol vector may be derived by (1) demodulating 
the symbol vector ?(«), possibly decoding and re-coding the demodulated data, and 
remodulating the demodulated data or re-coded data based on the signal constellations 
corresponding to the selected modulation schemes. The remodulated symbol vector, 



WO 03/050968 PCT/US02/38671 
020068 

31 

is then filtered by feed-back filter 518 with the frequency response of m 6 (/), 
and the output of filter 518 is provided to summer 516. 

|1099] Substituting equation (38) into equation (32) and assuming perfect decisions 
(i.e., !(w) = s(n)), the initial symbol estimate, ?(rc) , may be expressed as: 

s(Ai) = M / rs(/i) + M / |(«) , Eq(40) 
where |W = [$>(" + *,) + -1)L $(n)J . 

[1100] To determine the SNR associated with the initial symbol estimate, Jin), 
from the decision feedback equalizer, an unbiased minimum mean square error estimate 
is initially derived (similar to the MMSE-LE described above) by finding the 
conditional mean value of the transmitted symbol vector: 

£[!(*) I §(")] » M / fs(n) = G ^ s(n) , Eq(41) 

where G dfc = M f r = 9^ T . Next, the mean value of the z-th element of 
?(«), 7j(n) , is expressed as: 

E[s i (n)\s i (n)]=g dftii s i (n) 

where g dfe j7 is the i-th diagonal element of . 

[1101] To form the unbiased symbol estimate, s(n), similar to that described above 
for the MMSE-LE, a diagonal matrix whose elements are the inverse of the diagonal 
elements of G dfe is first defined as: 

Dcdfc =diag(l/g dfc>u , l/g dfe , 22 , A , \lg, Wr ) . Eq (42) 

The unbiased estimate, s(w) , can then be expressed as: 

§<») =DcdfelW = Dcdrc(M rW + Jljw) 

,~ v Eq(«) 

= Dc d rcM (rpCn) + I(n)) 



The resulting error covariance matrix is given by: 
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? M = W dfe =£{[s(/i)-D- , dfe s(/ I )] [s>)-s^/0D- dfe ] } 



Eq (44) 



= I-D-' dfc G 




[11 02) The SNR associated with the unbiased estimate, of the symbol 

transmitted on the z-th transmit antenna can then be expressed as: 



[1103] FIG. 5B is a block diagram of an embodiment of a decision feedback 
equalizer 322b, which is another embodiment of equalizer 322 in FIG. 3. Within 
decision feedback equalizer 322b, the received symbol vector, f (n) , from RX MIMO 
processor 160 is filtered by a feed-forward filter 534, which may implement the MMSE 
technique described above or some other linear spatial equalization technique. A 
summer 536 then combines the output from feed-forward filter 534 with the estimated 
distortion components from a feed-back filter 538 to provide the biased symbol 
estimate, s(n) , having the distortion component approximately removed. Initially, the 
estimated distortion components are zero and the symbol estimate, , is simply the 
output from filter 534. The initial estimate, sf(/i), for summer 536 is then multiplied 
with the matrix D~ dfc by a multiplier 540 to provide the unbiased estimate, of the 
transmitted symbol vector, s(n) . The unbiased estimate, s(n) , comprises the recovered 
symbol vector that is provided to RX data processor 162. 

[1104] Within RX data processor 162, symbol unmapping element 332 (in Fig. 3) 
provides demodulated data for the recovered symbol vector, s(n) . The demodulated 
data is then provided to a symbol mapping element 216x within DFE 322b and 
modulated to provide the remodulated symbol vector, l(n). Alternatively, the 
demodulated may be decoded, re-coded, and then provided to symbol mapping element 
216x. The remodulated symbols are estimates of the modulation symbols, s(«), 
transmitted from the transmitter. The remodulated symbol vector, is provided to 
feed-back filter 538, which filters the symbol vector to derive the estimated distortion 
components. Feed-back filter 538 may implement a linear spatial equalizer (e.g., a 
linear transversal equalizer). 



SNR. = 



1 



Eq(45) 



1 Sdfc.ii 
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[1105] For the DFE technique, the remodulated symbols are used to derive an 
estimate of the distortion generated by the already detected symbols. If the remodulated 
symbols are derived without errors (or with minimal errors), then the distortion 
component may be accurately estimated and the inter-symbol interference contributed 
by the already detected symbols may be effectively canceled out. The processing 
performed by feed-forward filter 534 and feed-back filter 538 is typically adjusted 
simultaneously to minimize the mean square error (MSE) of the inter-symbol 
interference in the recovered symbols. 

[1106] The DFE and MMSE techniques are described in further detail by S.L. 
Ariyavistakul et al in a paper entitled "Optimum Space-Time Processors with 
Dispersive Interference: Unified Analysis and Required Filter Span,'* IEEE Trans, on 
Communication, Vol. 7, No. 7, July 1999, and incorporated herein by reference. 

Maximum Likelihood Sequence Estimation 

[1107] Maximum likelihood sequence estimation (MLSE) for channels with inter- 
symbol interference (ISI) is performed by forming a set of path metrics for use in the 
Viterbi algorithm, which searches for the most likely transmitted sequence, given the 
observed received signal. The MLSE is described in further detail by Andrew J. Viterbi 
and Jim K. Omura in "Principles of Digital Communication and Coding/' McGraw-Hill, 
1979, which is incorporated herein by reference. 

[1108] The use of MLSE on wideband MIMO channels that have not been 
orthogonalized through eigenmode decomposition, however, is impractical due to the 
excessively high dimensionality of the channel state space. A Viterbi equalizer 
performing maximum likelihood sequence estimation for a MIMO channel has M rL 
states, where M is the size of the symbol alphabet, r<^N r is the number of independent 
transmitted data streams, and L is the channel memory. For example, in the simple case 
where QPSK signaling is used (M = 4) with four independent data streams (r = 4) and 
the channel has memory of one symbol (L = 1), the Viterbi equalizer would have 2 8 
states (i.e., 4 4 ] =2 8 ). 

[1109] Using time-domain eigenmode decomposition in conjunction with Viterbi 
MLSE greatly reduces the state space of the Viterbi equalizer. In this case, the received 
symbol streams can be equalized independently, so the size of the state space is now 
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linear in the number of independent data streams r, i.e., rM l . For the previous 
example, the state space would be reduced to 2 4 (i.e., 4 • 4 1 = 2 4 ). 
[1110] The objective of the MLSE approach is to choose the transmitted sequence 
of symbol vectors, s m (w) , that maximizes the metric: 

K m =Z{2Re[^(/i)O w W]-s;(/i)r w ^ of (/i)} . Eq(46) 

Since the blocks £(1 ) that make up T are diagonal, K m can be expressed as the sum of 
r metrics, each associated with one of the MEMO channel's time-domain eigen-modes: 

where 

*„('•) = 2») . Eq(48) 

n 

and 



,/» -2Keh(»)^s; if («-^r;(X)J 



Eq(49) 



L L 



k=0 ;=0 



[1111] The sequence metrics, k* w (0> ^ identical in form to sequence metrics 
associated with the MLSE for SISO channels with inter-symbol interference. 
Therefore, MLSE Viterbi equalization as is known in the art can be applied to the 
equalization of the individual received symbol streams, as follows. 
[1112] The path metric for received symbol stream i at stage n in the Viterbi 
algorithm is given by 

M l ( W ) = //:(/i) + M / (/i-l) . Eq(50) 



When sample n is received, (1) M values of ^{n) associated with each possible 
transmitted symbol, ^(n), are computed for each of the M L states associated with 
symbol steam / at sample time /i-l, and (2) M values of A/,.(/i), one associated with 
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each possible value of s mi (n) , are computed for each state. The largest value of Af,.(«) 

is then selected for each state at sample time n, and the sequence associated with this 
largest value is selected as the surviving sequence at that state. 

[1 113] Sequence decisions may be declared when path merge events occur, which is 
when all of the surviving sequences merge at a common previous state. Alternatively, 
sequence decisions may be declared when path truncation occur at a fixed delay, which 
may be used to force a choice when merge events have not yet occurred. 
[1 1 14] Several different types of equalizer have been described above, including the 
MMSE-LE, DFE, and MLSE. Each of these equalizers may be used to equalize the 
received symbols to provide the recovered symbols, which are estimates of the 
transmitted symbols. Other types of equalizer may also be used, and this is within the 
scope of the invention. By orthogonal izing the received symbol streams via eigenmode 
decomposition, the received symbol streams may be equalized independently, which 
can (1) greatly reduce the complexity of the equalizer selected for use and/or (2) allow 
for use of other types of equalizer that may be impractical otherwise. 
[1115] The techniques to transmit and receive data described herein may be 
implemented in various wireless communication systems, including but not limited to 
MMO and CDMA systems. These techniques may also be used for the forward link 
and/or the reverse link. 

[1116] The techniques described herein to process a data transmission at the 
transmitter and receiver may be implemented by various means. For example, these 
techniques may be implemented in hardware, software, or a combination thereof. For a 
hardware implementation, the elements used to perform various signal processing steps 
at the transmitter (e.g., to code and modulate the data, to derive the transmitter pulse- 
shaping matrices, to precondition the modulation symbols, and so on) or at the receiver 
(e.g., to derive the receiver pulse-shaping matrices, to precondition the received 
samples, to equalize the received symbols, to demodulate and decode the recovered 
symbols, and so on) may be implemented within one or more application specific 
integrated circuits (ASICs), digital signal processors (DSPs), digital signal processing 
devices (DSPDs), programmable logic devices (PLDs), field programmable gate arrays 
(FPGAs), processors, controllers, micro-controllers, microprocessors, other electronic 
units designed to perform the functions described herein, or a combination thereof. 
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[1117] For a software implementation, some or all of the signal processing steps at 
each of the transmitter and receiver may be implemented with modules (e.g., 
procedures, functions, and so on) that perform the functions described herein. The 
software codes may be stored in a memory unit (e.g., memories 132 and 172 in FIG. 1) 
and executed by a processor (e.g., controllers 130 and 170). The memory unit may be 
implemented within the processor or external to the processor, in which case it can be 
communicatively coupled to the processor via various means as is known in the art. 
[1118] Headings are included herein for reference and to aid in locating certain 
sections. These headings are not intended to limit the scope of the concepts described 
therein under, and these concepts may have applicability in other sections throughout 
the entire specification. 

[1119] The previous description of the disclosed embodiments is provided to enable 
any person skilled in the art to make or use the present invention. Various 
modifications to these embodiments will be readily apparent to those skilled in the art, 
and the generic principles defined herein may be applied to other embodiments without 
departing from the spirit or scope of the invention. Thus, the present invention is not 
intended to be limited to the embodiments shown herein but is to be accorded the widest 
scope consistent with the principles and novel features disclosed herein. 



[1 1 20] WHAT IS CLAIMED IS: 
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CLAIMS 

1. A method for transmitting data in a multiple-input multiple-output 
2 (MIMO) communication system, comprising: 

coding data in accordance with one or more coding schemes to provide coded 
4 data for a plurality of transmission channels in a MIMO channel; 

modulating the coded data in accordance with one or more modulation schemes 
6 to provide a plurality of streams of modulation symbols; 

deriving a pulse-shaping matrix based in part on an estimated response of the 
8 MIMO channel; 

preconditioning the plurality of modulation symbol streams based on the pulse- 
10 shaping matrix to derive a plurality of preconditioned signals; and 

transmitting the plurality of preconditioned signals over the MIMO channel. 

2 . The method of claim 1 , further comprising: 

2 determining an estimated channel response matrix for the MIMO channel; and 

decomposing the estimated channel response matrix to obtain a first sequence of 

4 matrices of eigen-vectors and a second sequence of matrices of singular values, and 

wherein the pulse-shaping matrix is derived based on the first and second 

6 sequences of matrices. 

3. The method of claim 2, wherein the estimated channel response matrix is 
2 given in the frequency domain and is decomposed in the frequency domain. 

4. The method of claim 2, wherein the estimated channel response matrix is 
2 decomposed based on singular value decomposition. 

5. The method of claim 2, wherein the estimated channel response matrix 
2 comprises a plurality of eigen-modes, and wherein eigen-modes associated with 

singular values below a particular threshold are not selected for use for data 
4 transmission. 
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6. The method of claim 2, wherein the singular values in each matrix in the 
2 second sequence are randomly ordered such that eigen-modes of the estimated channel 

response matrix are associated with approximately equal transmit power. 

7. The method of claim 2, further comprising: 

2 deriving a third sequence of matrices, based on the second sequence of matrices 

of singular values, having values indicative of transmit power allocated to eigen-modes 
4 of the estimated channel response matrix, and 

wherein the pulse-shaping matrix is derived based on the first and third 
6 sequences of matrices. 

8. The method of claim 7, wherein the third sequence of matrices is derived 
2 based on water-pouring analysis. 

9. The method of claim 1, wherein the pulse-shaping matrix comprises a 
2 plurality of sequences of time-domain values, and wherein the preconditioning is 

performed in the time domain by convolving the plurality of streams of modulation 
4 symbols with the pulse-shaping matrix. 

10. The method of claim 1, wherein the pulse-shaping matrix comprises a 
2 plurality of sequences of frequency-domain values, and wherein the preconditioning is 

performed in the frequency domain by multiplying a plurality of streams of transformed 
4 modulation symbols with the pulse-shaping matrix. 

11. The method of claim 1, wherein the pulse-shaping matrix is derived to 
2 maximum capacity by allocating more transmit power to transmission channels with 

higher signal-to-noise-and-interference ratios (SNRs). 

12. The method of claim 1, wherein the pulse- shaping matrix is derived to 
2 provide approximately equal received signal-to-noise-and-interference ratios (SNRs) for 

the plurality of modulation symbol streams. 



WO 03/050968 PCT/US02/38671 
020068 

39 

13. The method of claim 1, wherein separate coding and modulation 
schemes are used for each transmission channel. 



14. The method of claim 1, wherein a common coding and modulation 
2 scheme is used for all transmission channels. 



15. A method for transmitting data in a multiple-input multiple-output 
2 (MIMO) communication system, comprising: 

coding data in accordance with one or more coding schemes to provide coded 
4 data for a plurality of transmission channels in a MIMO channel; 

modulating the coded data in accordance with one or more modulation schemes 
6 to provide a plurality of streams of modulation symbols; 

determining an estimated channel response matrix for the MEMO channel; 
8 decomposing the estimated channel response matrix to obtain a first sequence of 

matrices of eigen-vectors and a second sequence of matrices of singular values; 
10 deriving a third sequence of matrices, based on the second sequence of matrices 

of singular values, having values indicative of transmit power allocated to eigen-modes 
1 2 of the estimated channel response matrix; 

deriving a pulse-shaping matrix based on the first and third sequences of 
14 matrices; 

preconditioning the plurality of modulation symbol streams based on the pulse- 
16 shaping matrix to derive a plurality of preconditioned signals; and 

transmitting the plurality of preconditioned signals over the MIMO channel. 

16. A memory communicatively coupled to a digital signal processing 
2 device (DSPD) capable of interpreting digital information to: 

code data in accordance with one or more coding schemes to provide coded data 
4 for a plurality of transmission channels in a MIMO channel; 

modulate the coded data in accordance with one or more modulation schemes to 
6 provide a plurality of streams of modulation symbols; 

derive a pulse-shaping matrix based in part on an estimated response of the 
8 MIMO channel ; and 
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precondition the plurality of modulation symbol streams based on the pulse- 
10 shaping matrix to derive a plurality of preconditioned signals for transmission over the 
MIMO channel 

17. A method for receiving a data transmission in a multiple-input multiple- 
2 output (MIMO) communication system, comprising: 

determining an estimated channel response matrix for a MIMO channel used for 
4 the data transmission; 

decomposing the estimated channel response matrix to obtain a first sequence of 
6 matrices of eigen-vectors; 

deriving a pulse-shaping matrix based on the first sequence of matrices; and 
8 preconditioning a plurality of received signals based on the pulse-shaping matrix 

to obtain a plurality of streams of received symbols. 

18. The method of claim 17, wherein the preconditioning is performed in the 
2 time domain based on a time-domain pulse-shaping matrix. 

19. The method of claim 17, wherein the preconditioning is performed in the 
2 frequency domain and includes 

transforming the plurality of received signals to the frequency domain; 

4 multiplying the transformed received signals with a frequency-domain pulse- 

shaping matrix to derived a plurality of preconditioned signals; and 

6 transforming the plurality of preconditioned signals to the time domain to obtain 

the plurality of received symbol streams. 

20. The method of claim 17, wherein the preconditioning orthogonalizes the 
2 plurality of streams of received symbols. 



2 



21 . The method of claim 1 7, further comprising: 

equalizing the plurality of received symbol streams to derive a plurality of 
recovered symbol streams. 
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22. The method of claim 21, wherein the equalizing is performed separately 
2 for each received symbol stream. 



23. The method of claim 21, wherein the equalizing is performed based on a 
2 minimum mean square error linear equalizer (MMSE-LE). 

24. The method of claim 21, wherein the equalizing is performed based on a 
2 decision feedback equalizer (DFE). 

25. The method of claim 21 , wherein the equalizing is performed based on a 
2 maximum likelihood sequence estimation (MLSE) equalizer. 

26. The method of claim 21, further comprising: 

2 demodulating the plurality of recovered symbol streams in accordance with one 

or more demodulation schemes to provide a plurality of demodulated data streams; and 

4 decoding the plurality of demodulated data streams in accordance with one or 

more decoding schemes to provide decoded data. 

27. The method of claim 17, further comprising: 

2 deriving channel state information (CSI) comprised of the estimated channel 

response matrix and signal-to-noise-and-interference ratios (SNRs) for a plurality of 
4 transmission channels of the MEMO channel; and 

sending the CSI back to a transmitter of the data transmission. 

28. A method for receiving a data transmission in a multiple-input multiple- 
2 output (MIMO) communication system, comprising: 

determining an estimated channel response matrix for a MIMO channel used for 
4 the data transmission; 

decomposing the estimated channel response matrix to obtain a first sequence of 
6 matrices of ei gen- vectors; 

deriving a pulse-shaping matrix based on the first sequence of matrices; 
8 preconditioning a plurality of received signals based on the pulse-shaping matrix 

to obtain a plurality of streams of received symbols; 
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10 equalizing the plurality of received symbol streams to derive a plurality of 

recovered symbol streams; 
12 demodulating the plurality of recovered symbol streams in accordance with one 

or more demodulation schemes to provide a plurality of demodulated data streams; and 
14 decoding the plurality of demodulated data streams in accordance with one or 

more decoding schemes to provide decoded data. 

29. A memory communicatively coupled to a digital signal processing 
2 device (DSPD) capable of interpreting digital information to: 

determine an estimated channel response matrix for a MIMO channel used for 
4 the data transmission; 

decompose the estimated channel response matrix to obtain a first sequence of 
6 matrices of eigen-vectors; 

derive a pulse-shaping matrix based on the first sequence of matrices; and 
8 precondition a plurality of received signals based on the pulse-shaping matrix to 

obtain a plurality of streams of received symbols. 

30. A transmitter unit in a multiple-input multiple-output (MIMO) 
2 communication system, comprising: 

a TX data processor operative to code data in accordance with one or more 
4 coding schemes to provide coded data for a plurality of transmission channels in a 
MIMO channel, and to modulate the coded data in accordance with one or more 
6 modulation schemes to provide a plurality of streams of modulation symbols; 

a TX MIMO processor operative to derive a pulse-shaping matrix based in part 
8 on an estimated response of the MIMO channel, and to precondition the plurality of 
modulation symbol streams based on the pulse-shaping matrix to provide a plurality of 
10 preconditioned signals; and 

one or more transmitters operative to condition and transmit the plurality of 
1 2 preconditioned signals over the MIMO channel. 

31. The transmitter unit of claim 30, wherein the TX MIMO processor is 
2 further operative to determine an estimated channel response matrix for the MIMO 

channel, decompose the estimated channel response matrix to obtain a first sequence of 
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4 matrices of eigen-vectors and a second sequence of matrices of singular values, and 
derive the pulse-shaping matrix based on the first and second sequences of matrices. 

32. The transmitter unit of claim 31, wherein the TX MIMO processor is 
2 further operative to decompose the estimated channel response matrix in the frequency 
domain using singular value decomposition. 



33. The transmitter unit of claim 31, wherein the TX MIMO processor is 
2 further operative to derive, based on the second sequence of matrices of singular values, 

a third sequence of matrices of values indicative of transmit power allocated to eigen- 
4 modes of the estimated channel response matrix, and to derive the pulse-shaping matrix 
based on the first and third sequences of matrices. 

34. A transmitter apparatus in a multiple-input multiple-output (MIMO) 
2 communication system, comprising: 

means for coding data in accordance with one or more coding schemes to 
4 provide coded data for a plurality of transmission channels in a MIMO channel; 

means for modulating the coded data in accordance with one or more 
6 modulation schemes to provide a plurality of streams of modulation symbols; 

means for deriving a pulse-shaping matrix based in part on an estimated 
8 response of the MIMO channel; 

means for preconditioning the plurality of modulation symbol streams based on 
1 0 the pulse-shaping matrix to derive a plurality of preconditioned signals; and 

means for transmitting the plurality of preconditioned signals over the MIMO 
12 channel. 

35. A digital signal processor for use in a multiple-input multiple-output 
2 (MIMO) communication system, comprising: 

means for coding data in accordance with one or more coding schemes to 
4 provide coded data for a plurality of transmission channels in a MIMO channel; 

means for modulating the coded data in accordance with one or more 
6 modulation schemes to provide a plurality of streams of modulation symbols; 
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means for deriving a pulse-shaping matrix based in part on an estimated 
8 response of the MIMO channel; and 

means for preconditioning the plurality of modulation symbol streams based on 
1 0 the pulse-shaping matrix to derive a plurality of preconditioned signals. 



36. A receiver unit in a multiple-input multiple-output (MIMO) 
2 communication system, comprising: 

a RX MIMO processor operative to determine an estimated channel response 
4 matrix for a MIMO channel used for data transmission, decompose the estimated 
channel response matrix to obtain a first sequence of matrices of eigen-vectors, derive a 
6 pulse-shaping matrix based on the first sequence of matrices, and precondition a 
plurality of received signals based on the pulse-shaping matrix to obtain a plurality of 
8 streams of received symbols; and 

a RX data processor operative to demodulate the plurality of received symbol 
10 streams in accordance with one or more demodulation schemes, and to decode a 
plurality of demodulated data streams in accordance with one or more decoding 
1 2 schemes to provide decoded data. 

37. The receiver unit of claim 36, wherein the RX MIMO processor includes 
2 an equalizer operative to equalize the plurality of received symbol streams to 

provide a plurality of recovered symbol streams, and 
4 wherein the RX data processor is operative to demodulate and decode the 

plurality of recovered symbol streams to provide the decoded data. 



38. The receiver unit of claim 37, wherein the equalizer is a minimum mean 
2 square error linear equalizer (MMSE-LE). 

39. The receiver unit of claim 37, wherein the equalizer is a decision 
2 feedback equalizer (DFE). 



40. The receiver unit of claim 37, wherein the equalizer is a maximum 
2 likelihood sequence estimation (MLSE) equalizer. 
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41. The receiver unit of claim 37, wherein the equalizer is operative to 
separately equalize each received symbol stream. 



42. A receiver apparatus in a multiple-input multiple-output (MIMO) 
2 communication system, comprising: 

means for determining an estimated channel response matrix for a MIMO 
4 channel used for the data transmission; 

means for decomposing the estimated channel response matrix to obtain a first 
6 sequence of matrices of eigen- vectors; 

means for deriving a pulse-shaping matrix based on the first sequence of 
8 matrices; and 

means for preconditioning a plurality of received signals based on the pulse- 
10 shaping matrix to obtain a plurality of streams of received symbols. 

43. A digital signal processor in a multiple-input multiple-output (MIMO) 
2 communication system, comprising: 

means for determining an estimated channel response matrix for a MIMO 
4 channel used for the data transmission; 

means for decomposing the estimated channel response matrix to obtain a first 
6 sequence of matrices of eigen-vectors; 

means for deriving a pulse-shaping matrix based on the first sequence of 
8 matrices; and 

means for preconditioning a plurality of received signals based on the pulse- 
10 shaping matrix to obtain a plurality of streams of received symbols. 
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